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The SciBooNE Collaboration has performed a search for charged current coherent pion production 
from muon neutrinos scattering on carbon, u^^^C — > /^"^^CTr'*', with two distinct data samples. No 
evidence for coherent pion production is observed. We set 90% confidence level upper limits on the 
cross section ratio of charged current coherent pion production to the total charged current cross 
section at 0.67 x lO^'^ at mean neutrino energy 1.1 GeV and 1.36 x lO"'^ at mean neutrino energy 
2.2 GeV. 

PACS numbers: 13.15.+g, 13.60.Le, 25.30.Pt, 95.55.Vj 



I. INTRODUCTION 



Although they have been studied for decades, neutrino- 
nucleus cross sections between 100 MeV and 10 GcV 
energy are still known with very poor accuracy. The 
demand for precise cross section measurements in this 
energy regime is driven by the needs of the next genera- 
tion of neutrino oscillation experiments in their pursuit of 
sub-leading flavor oscillation and charge-parity violation 
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[ij, l^l . Several interaction channels contribute to the to- 
tal neutrino-nucleus cross section in the neutrino energy 
range of a few GeV. Interactions producing single pions 
(charged or neutral) account for a large cross section frac- 
tion, which must be understood because they form sig- 
nificant backgrounds for neutrino oscillation searches. 

It has been known for years that neutrinos can produce 
pions by interacting coherently with the nucleons forming 
the target nucleus. The cross section for this process is 
expected to be smaller than incoherent pion production, 
the latter being dominated by neutrino-induced baryonic 
resonance excitation off a single nucleon bound in a nu- 
cleus. Moreover, coherent pion production is compara- 
tively poorly understood, although it is characterized by 
a distinct signature consisting of a nucleus left in the 
ground state (no nuclear breakup occurs) and a forward 
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scattered pion. Both charged current and neutral cur- 
rent coherent modes are possible, v^A — > fj.^ Att^ and 
v^A ^ i/^Att'^, where A is a nucleus. 

Several theoretical models describing coherent pion 
production have been proposed, using different for- 
malisms to describe the relevant physics. A first class 
of models is built on the basis of Adler's PCAC the- 
orem 0], relating the neutrino- nucleus cross section to 
that of a pion interacting with a nucleus at = 0; the 
extrapolation to Q'^ ^ is performed via a propagator 
term 0, H, @, 0, S]- A second commonly- used formal- 
ism is based on the description of the coherent produc- 
tion of A resonances on nuclei by using a modified A- 
propagator and a distorted wave-function for the pion 
@j [3j [m, [13 • While the relationship between neutral 
current and charged current modes, and that between 
neutrino and antineutrino coherent pion production cross 
sections, are relatively well known, order-of-magnitude 
variations on absolute coherent pion production cross sec- 
tions are expected within these models. In addition, the 
cross section dependence on neutrino energy and on tar- 
get material is also uncertain. It is therefore imperative 
that more experimental input on coherent pion produc- 
tion in neutrino-nucleus interactions is gathered in the 
near future. 

Coherent pion production in neutrino-nucleus interac- 
tions has already been the subject of several experimen- 
tal campaigns. The neutrino energy range between 1 
and 100 GeV has been investigated, including both the 
charged current and neutral current modes, and using 
both neutrino and antineutrino probes. A result that has 
drawn much attention in the neutrino physics community 
has been the recent non-observation of charged current 
coherent pion production by the K2K experiment with 
a 1.3 GeV wide-band neutrino beam [l3|. This is moti- 
vated by the fact that the K2K Collaboration has quoted 
an upper limit for the ratio of the charged current co- 
herent pion production cross section to the charged cur- 
rent inclusive cross section that is well below the predic- 
tion of the original Rein-Sehgal model Q that has been 
adopted in the past to describe coherent pion production 
processes. In addition, even within more recently pro- 
posed models, it is often difficult to reconcile this new 
and accurate null result at low energies with previous 
measurements. On the one hand, evidence for neutral 
current coherent pion production in a neutrino energy 
range that is similar to K2K has been unambiguously 
reported from the Aachen-Padova [l^ and Gargamelle 
experimental data [ist first, and more recently also by 
the MiniBooNE Collaboration On the other hand, 
while no measurements of charged current coherent pion 
production other than the K2K one exist in a similar neu- 
trino energy range, there exist charged current coherent 
pion production positive results at higher energies (7- 
100 GeV neutrino energy) from the SKAT [13], CHARM 
[3, BEBC [H,!!!, and FNAL E632 ^ experiments. 

In this paper, we discuss the first measurement of 
charged current coherent pion production by neutrinos 



in the SciBooNE experiment This is a particularly 
interesting test of the K2K null result, probing a similar 
neutrino energy range and same target material. Also, 
compared to K2K, SciBooNE 's result presented here is 
based on a higher-statistics data sample and uses an im- 
proved analysis, as will be described below. 

This paper is organized as follows. Section [Til describes 
the neutrino beam-line and the neutrino flux simulation. 
The simulation of neutrino interactions with nuclei are 
described in Section Hill The detector configuration and 
simulation are described in Section IIVI A summary of 
the data set and experimental performance is given in 
Section |Vl The data analysis, including the event selec- 
tion and Monte Carlo (MC) tuning, is described in detail 
in Section rvTl The results of the analysis and discussion 
are presented along with a summary of systematic un- 
certainties in Section IVIIi and the final conclusions are 
given in Section IVlIII 



II. NEUTRINO BEAM 

The SciBooNE detector has been exposed to the 
Booster Neutrino Beam (BNB) located at Fermilab. The 
BNB is a high-intensity, conventional neutrino beam 
which has been serving the MiniBooNE experiment since 
2002. 



A. Beam-line Description 

The primary beam uses protons accelerated to 8 GeV 
kinetic energy by the Fermilab Booster. Selected batches 
containing approximately 4-5 x 10^^ protons are extracted 
and bent toward the BNB target hall via dipole magnets. 
Each spill is composed of 81 bunches of protons, approx- 
imately 6 ns wide each and 19 ns apart, for a total spill 
duration of 1.6 ^s. 

Beam proton trajectories and positions are monitored 
on a pulse-by-pulse basis. The typical beam alignment 
and divergence measured by the beam position monitors 
located near the target are within 1 mm and 1 mrad of the 
nominal target center and axis direction, respectively; the 
typical beam focusing on target measured by beam profile 
monitors is of the order of 1-2 mm (RMS) in both the 
horizontal and vertical directions. These parameters are 
well within the experiment requirements. The number of 
protons delivered to the BNB target is measured for each 
proton batch using two toroids located near the target 
along the beam-line. The toroid calibration, performed 
on a pulse-by-pulse basis, provides a measurement of the 
number of protons to BNB with a 2% accuracy. 

Primary protons from the 8 GeV beamline strike a 
thick beryllium target located in the BNB target hall. 
Hadronic interactions of the protons with the target ma- 
terial produce a beam of secondary mesons (pions and 
kaons). The target is made of seven cylindrical slugs for 
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a total target length of 71.1 cm, or about 1.7 inelastic 
interaction lengths. 

The beryllium target is surrounded by a magnetic fo- 
cusing horn, bending and sign-selecting the secondary 
particles that emerge from the interactions in the tar- 
get along the direction pointing to the SciBooNE detec- 
tor. The focusing is produced by the toroidal magnetic 
field present in the air volume between the horn's two 
coaxial conductors made of aluminum alloy. The horn 
current pulse is approximately a half-sinusoid of ampli- 
tude 174 kA, 143 yus long, synchronized to each beam 
spill. The polarity of the horn current flow can be (and 
has been) switched, in order to focus negatively-charged 
mesons, and therefore produce an antineutrino instead of 
a neutrino beam. 

The beam of focused, secondary mesons emerging from 
the target/horn region is further collimated via passive 
shielding, and allowed to decay into neutrinos in a cylin- 
drical decay region filled with air at atmospheric pressure, 
50 m long and 90 cm in radius. A beam absorber located 
at the end of the decay region stops the hadronic and 
muonic component of the beam, and only a pure neu- 
trino beam pointing toward the detector remains, mostly 
from 7r+ — » /x+i/^ decays. 



B. Neutrino Flux Prediction 

Neutrino flux predictions at the SciBooNE detector 
location are obtained via a GEANT4-based [1^ beam 
Monte Carlo simulation. The same simulation code de- 
veloped by the MiniBooNE Collaboration is used [23 |. 

In the simulation code, a realistic description of the 
geometry and materials present in the BNB target hall 
and decay region is used. Primary protons are gener- 
ated according to the expected beam optics properties 
upstream of the target. The interactions of primary pro- 
tons with the beryllium target are simulated according 
to statc-of-thc-art hadron interaction data. Of particular 
importance for this analysis is 7r+ production in proton- 
beryllium interactions, which uses experimental input 
from the HARP [25i] and BNL E910 ^ experiments. 
Production of secondary protons, neutrons, charged pi- 
ous, and charged and neutral kaons is taken into account, 
and elastic and quasi-elastic scattering of protons in the 
target are also simulated. Particles emanating from the 
primary proton-beryllium interaction in the target are 
then propagated within the GEANT4 framework, which 
accounts for all relevant physics processes. Hadronic re- 
interactions of pions and nucleons with beryllium and 
aluminum materials are particularly important and are 
described by custom models, while other hadronic pro- 
cesses and all electromagnetic processes (energy loss, 
multiple scattering, effect of horn magnetic field, etc.) 
are described according to default GEANT4 physics lists. 
A second, FORTRAN-based Monte Carlo code uses the 
output of the GEANT4 program as input, and is respon- 
sible for generating the neutrino kinematics distributions 



from meson and muon decays, and for obtaining the fi- 
nal neutrino fluxes extrapolated to the SciBooNE detec- 
tor with negligible beam Monte Carlo statistical errors. 
Current best knowledge of neutrino-producing meson and 
muon decay branching fractions, and decay form factors 
in three-body semi-leptonic decays, are used. Polariza- 
tion effects in muon decays are also accounted for. 

Once produced by the simulation, neutrinos are ex- 
trapolated along straight lines toward the SciBooNE de- 
tector. All neutrinos whose ray traces cross any part of 
the detector volume are considered for SciBooNE flux 
predictions. Based on accurate survey data, the distance 
between the center of the beryllium target and the cen- 
ter of the SciBar detector is taken to be 99.9 m, with the 
SciBooNE detector located on beam axis within a toler- 
ance of a few cm. Each simulated neutrino interaction 
is linked to its detailed beam information and history, 
which includes neutrino flavor, energy, parent type and 
kinematics, and ray trace entry and exit points within 
the detector volume; the ray trace information is used to 
determine the incoming neutrino's direction and inter- 
action location. Proper weights for each beam neutrino 
event are computed, using this beam neutrino informa- 
tion, as well as information from the interaction and de- 
tector simulation: neutrino interaction probability, and 
detailed SciBooNE detector geometry and specifications. 




Ev (GeV) 



FIG. 1: Neutrino flux prediction at the SciBooNE detector as 
a function of neutrino energy normalized per unit area, 
proton on target (POT) and neutrino energy bin width. The 
spectrum is averaged within 2.12 m from the beam center. 
The total flux and contributions from individual neutrino fla- 
vors are shown. 

The neutrino flux prediction at the SciBooNE detector 
location and as a function of neutrino energy is shown 
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in Fig. [TJ A total neutrino flux per proton on target 
of 2.2 X 10^® cm"^ is expected at the SciBooNE detec- 
tor location and in neutrino running mode (positive horn 
polarity), with a mean neutrino energy of 0.7 GeV. The 
flux is dominated by muon neutrinos (93% of total), with 
small contributions from muon antincutrinos (6.4%), and 
electron neutrinos and antincutrinos (0.6% in total). For 
the neutrino flux predictions used in this analysis, no 
information from BNB (SciBooNE or MiniBooNE) neu- 
trino data is used as experimental input. 

III. NEUTRINO INTERACTION SIMULATION 

The neutrino interactions with nuclear targets are sim- 
ulated with the NEUT program library [23, which 
is used in the Kamiokande, Super-Kamiokande, K2K, 
and T2K experiments. NEUT handles protons, oxy- 
gen, carbon, and iron as nuclear targets in the energy 
range from 100 McV to 100 TcV. In NEUT, the follow- 
ing neutrino interactions in both neutral and charged 
currents are simulated: quasi-elastic scattering (vN — > 
(N'), single meson production {lyN £N'm), single 
gamma production (i^N — > iN'j), coherent tt produc- 
tion (i^i2c(or^^Fe) h i2c(or^'^Fe)), and deep inelastic 
scattering {vN iN'hadrons), where A'' and N' are the 
nucleons (proton or neutron), £ is the lepton, and m is 
the meson. Following the primary neutrino interactions 
in nuclei, re-interactions of the mesons and hadrons with 
the nuclear medium are also simulated. 



A. Coherent tt production 

The signal for this analysis, coherent pion production, 
is a neutrino interaction with a nucleus which remains 
intact, releasing one pion with the same charge as the 
incoming weak current. Because of the small momen- 
tum transfer to the target nucleus the outgoing lepton 
and pion tend to go in the forward direction (in the lab 
frame) . The formalism developed by Rein and Sehgal ^ 
is used to simulate the interactions, including the recent 
correction of lepton mass effects [g. The axial vector 
mass. Ma, is set to 1.0 GeV/c^. The nuclear radius pa- 
rameter i?o is set to 1.0 fm. For the total and inelastic 
pion-nucleon cross sections used in the formalism, the 
fitted results given in Rein and Schgal's paper are em- 
ployed. The total cross section on ^^C is shown in Fig. [2l 
with comparisons of other models discussed in the intro- 
duction. The Rein and Sehgal model predicts charged 
current coherent pion production to be approximately 
1% of the total neutrino interactions in SciBooNE. 



B. Quasi-elastic scattering 

The dominant interaction in the SciBooNE neutrino 
energy range is quasi-elastic scattering, which is imple- 
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FIG. 2: Cross section for i/^^^C— > /^"tt^^'^C interaction. The 
solid line represents the Rein and Sehgal model with lepton 
mass effects 0](the default model of the signal for this anal- 
ysis), the dashed line represents the Rein and Sehgal model 
without lepton mass effectsja] , the dotted line represents the 
model of Kartavtsev et al. and the dashed-dotted line rep- 
resents the model of Alvarez-Ruso et al. The model of 
Singh et al. gives a cross section similar to the model of 
Alvarez-Ruso et al. 

mcnted using the model of Llewellyn-Smith [2^. For 
scattering off nucleons in the nucleus, we use the rela- 
tivistic Fermi gas model of Smith and Moniz [sO]. The 
nucleons are treated as quasi-free particles and the Fermi 
motion of nucleons along with the Pauli exclusion prin- 
ciple is taken into account. The momentum distribution 
of the target nucleon is assumed to be flat up to a fixed 
Fermi surface momentum of 217 MeV/c for carbon and 
250 MeV/c for iron. The same Fermi momentum distri- 
bution is also used for all of the other nuclear interactions. 
The nuclear potential is set to 27 MeV for carbon and 
32 MeV for iron. Both vector and axial-vector form fac- 
tor are assumed to be dipole. The vector mass in quasi- 
elastic scattering is set to be 0.84 GeV/c^. The axial 
vector mass. Ma, is set to be 1.21 GeV/c^ as suggested 
by recent results [3l|, [13] . 



C. Single meson production via baryon resonances 

The second most probable interaction in SciBooNE is 
the resonant single meson production of tt, K, and rj de- 
scribed by the model of Rein and Sehgal [33| . The model 
assumes an intermediate baryon resonance, A''*, in the 
reaction of -> eN*,N* N'm. The differential 
cross section of single meson production depends on the 
amplitude for the production of a given resonance and 
the probability of the baryon resonance decay to the me- 
son. All intermediate baryon resonances with mass less 
than 2 GeV/c^ are included. Those baryon resonances 
with mass greater than 2 GeV/c^ are simulated as deep 
inelastic scattering. Lepton mass effects from the non- 
conservation of lepton current and the pion-pole term 
in the hadronic axial vector current are included in the 
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simulation [3J, |35|. 

To determine the angular distribution of a pion in the 
final state, Rein's method [IHl is used for the P33(1232) 
resonance. For other resonances, the directional distri- 
bution of the generated pion is set to be isotropic in the 
resonance rest frame. The angular distribution of 7r+ 
has been measured for v^p fi~pTT^ [s^l and the re- 
sults agree well with NEUT's prediction. Pauli blocking 
is accounted for in the decay of the baryon resonance by 
requiring the momentum of the nucleon to be larger than 
the Fermi surface momentum. Pion-less A decay is also 
taken into account, where 20% of the events do not have a 
pion and only the lepton and nucleon are generated [11] . 
The axial vector mass, Ma, is set to be 1.21 GeV/c^. 



D. Deep inelastic scattering 

The cross section for deep inelastic scattering (DIS) 
is calculated using the GRV98 parton distribution func- 
tions [SE^ . Additionally, we have included the corrections 
in the small region developed by Bodek and Yang [40l |. 
In the calculation, the hadronic invariant mass, W, is 
required to be larger than 1.3 GeV/c^. Also, the mul- 
tiplicity of pions is restricted to be larger than or equal 
to two for 1.3 < W < 2.0 GeV/c^, because single pion 
production is already included in the simulation, as de- 
scribed above. The multi-hadron final states are sim- 
ulated with two models: a custom-made program (4l| 
for the event with W between 1.3 and 2.0 GeV/c^ and 
PYTHIA/JETSET ^ for the events with W larger than 
2 GeV/c2. 



E. Intra-nuclear interactions 

The intra-nuclear interactions of mesons and nucleons 
produced in neutrino interactions in the nuclei are sim- 
ulated. These interactions are treated using a cascade 
model, and each of the particles is traced until it escapes 
from the nucleus. 

Among all the interactions of mesons and nucleons, the 
interactions of pions are most important to this analysis. 
The inelastic scattering, charge exchange and absorption 
of pions in the nuclei are simulated. The interaction cross 
sections of pions in the nuclei are calculated using the 
model by Salcedo et al. [i^, which agrees well with past 
experimental data [ij. If inelastic scattering or charge 
exchange occurs, the direction and momentum of pions 
are determined by using results from a phase shift anal- 
ysis of pion- nucleus scattering experiments [i^. When 
calculating the pion scattering amplitude, Pauli blocking 
is taken into account by requiring the nucleon momentum 
after the interaction to be larger than the Fermi surface 
momentum at the interaction point. 

Re-interactions of the recoil protons and neutrons pro- 
duced in neutrino interactions are also important, be- 
cause the proton tracks are used to classify the neutrino 



event type. Nucleon-nucleon interactions modify the out- 
going nucleon's momentum and direction. Both elastic 
scattering and pion production are considered. In or- 
der to simulate these interactions, a cascade model is 
again used and the generated particles in the nucleus are 
tracked using the same code as for the mesons. 

No de-excitation gamma-ray from the carbon nucleus 
is simulated when nuclear breakup occurs. 



IV. NEUTRINO DETECTOR 

The SciBooNE detector is located 100 m downstream 
from the beryllium target on the axis of the beam. The 
detector comprises three sub-detectors: a fully active and 
finely segmented scintillator tracker (SciBar), an electro- 
magnetic calorimeter (EC), and a muon range detector 
(MRD). 



A. Detector Description 

Fig. [3] shows an event display of a typical muon neu- 
trino charged current single charged pion event candi- 
date. Detector coordinates are shown in the figure. Sci- 
BooNE uses a right-handed Cartesian coordinate system 
in which the z axis is the beam direction and the y axis 
is the vertical upward direction. The origin is located 
on the most upstream surface of SciBar in the z dimen- 
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FIG. 3: Event display of a typical muon neutrino charged cur- 
rent single charged pion event candidate in SciBooNE data. 
The neutrino beam runs from left to right in this figure, en- 
countering SciBar, the EC and MRD, in that order. The 
circles on SciBar and the EC indicate ADC hits for which the 
area of the circle is proportional to the energy deposition in 
that channel. Filled boxes in the MRD show ADC hits in 
time with the beam window. 



6 



sion, and at the center of the SciBar scintillator plane 
in the x and y dimensions. Since each sub-detector is 
read out both vertically and horizontally, two views arc 
defined; the top view [z-x projection) and the side view 
{z-y projection). 

The SciBar detector [iHl is positioned upstream of the 
other sub-detectors. The primary role of SciBar is to 
reconstruct the neutrino-nucleus interaction vertex and 
detect charged particles produced by neutrino interac- 
tions. Moreover, SciBar is capable of particle identifica- 
tion based on deposited energy. SciBar was designed and 
built as a near detector for the K2K experiment. After 
K2K's completion, SciBar was relocated to the Fermilab 
BNB for SciBooNE. 

SciBar consists of 14,336 extruded plastic scintillator 
strips which serve as the target for the neutrino beam 
as well as the active detection medium. Originally pro- 
duced by Fermilab, each strip has a dimension of 1.3 x 
2.5 X 300 cm'^. The scintillators are arranged vertically 
and horizontally to construct a 3 x 3 x 1.7 m^ volume 
with a total mass of 15 tons. Each strip is read out by a 
wavelength shifting (WLS) fiber attached to a 64-channcl 
multi- anode PMT. Charge and timing information from 
each PMT is recorded by front-end electronics boards 
(FEB) attached directly to the PMT and a back-end 
VME module [13 . The FEB uses VA/TA ASICs; the 
VA handles charge information from the PMT with a 32- 
channcl preamplifier chip with a shaper and multiplexer, 
while the TA provides timing information by taking the 
logical "OR" of 32 channels. The charge and timing infor- 
mation are digitized by ADC and multi-hit TDC modules 
on back-end electronics, and read out through the VME 
bus. 

The gains of all PMT channels were measured prior 
to installation in K2K. SciBar is eqiuppcd with a gain 
calibration system comprised of LEDs to monitor and 
correct gain drift during the data taking; the gain stabil- 
ity is monitored with precision better than 1%. Cosmic 
ray data arc also employed to calibrate the PMT gains 
and scintillator light yield, including attenuation of the 
WLS-fibers. These calibration data, LED and cosmic, 
are taken between beam spills and continuously moni- 
tored. Calibration data verify the light yield was sta- 
ble within 1% during operation. The timing resolution 
for minimum-ionizing particles was evaluated with cos- 
mic ray data to be 1.6 ns. The average light yield for 
minimum-ionizing particles is approximately 20 photo- 
electrons per 1.3 cm path length, and the typical pedestal 
width is below 0.3 photoelectron. The hit finding ef- 
ficiency evaluated with cosmic ray data is more than 
99.8%. The minimum length of a reconstructable track 
is approximately 8 cm (three layers hit in each view). 
The track finding efficiency for single tracks of 10 cm or 
longer is more than 99%. 

The EC is located just downstream of SciBar, and is 
designed to measure the electron neutrino contamination 
in the beam and tag photons from 7r° decay. The EC is 
a "spaghetti" type calorimeter comprised of 1 mm di- 



ameter scintillating fibers embedded in lead foil. The 
calorimeter is made of modules of dimensions 262 x 8 x 
4 cm'^. Each module is read out by two green-extended 
1 inch Hamamatsu PMTs per side, 256 PMTs total. The 
modules were originally built for the CHORUS neutrino 
experiment at CERN [4^ and later used in HARP and 
then K2K. The modules construct one vertical and one 
horizontal plane, and each plane has 32 modules. The EC 
has a thickness of 11 radiation lengths along the beam di- 
rection. The planes cover an active area of 2.7 x 2.6 m^. 

The charge information from each PMT is recorded. 
A minimum ionizing particle with a minimal path 
length deposits approximately 91 MeV in the EC. The 
energy resolution for electrons was measured to be 
IA%/^JE (GeV) using a test beam [H. 

The MRD is installed downstream of the EC and is 
designed to measure the momentum of muons produced 
by charged-current neutrino interactions. The MRD was 
constructed for SciBooNE at Fermilab, primarily out of 
parts recycled from past experiments. It has 12 iron 
plates with thickness of 5 cm which are sandwiched be- 
tween planes of 6 mm thick scintillation counters, 13 al- 
ternating horizontal and vertical planes, which are read 
out via 2 inch PMTs from a variety of past experiments; 
there are 362 PMTs total. Each iron plate covers an area 
of 274 X 305 cm^ . The total mass of absorber material 
is approximately 48 tons. The MRD measures the mo- 
mentum of muons up to 1.2 GeV/c using the observed 
muon range. Charge and timing information from each 
PMT are recorded. Hit finding efficiency was continu- 
ously monitored using cosmic ray data taken between 
beam spills; the average hit finding efficiency is 99%. 

SciBooNE has two global triggers, the beam trigger 
and the off-beam trigger. Two types of data are collected 
in one beam cycle, neutrino data with the beam trigger 
and calibration data with the off-beam trigger. One cycle 
is about 2 sec which is defined by the accelerator timing 
sequence. The BNB receives one train of proton beam 
pulses per cycle, with a maximum of 10 pulses in a row 
at 15 Hz. 

A fast timing signal sent by the extraction magnet on 
BNB pulses establishes a beam-trigger. Once the beam 
trigger condition is set, all sub-detector systems read out 
all channels irrespective of hit occupancy (i.e. whether or 
not a neutrino interaction occurred), ensuring unbiased 
neutrino data. 

After the beam trigger turns off, the off-beam trig- 
ger condition is automatically set and each sub-detector 
takes calibration data. There are three types of calibra- 
tion data: pedestal, LED (only for SciBar) and cosmic 
ray data. The pedestal and LED data are collected once 
per cycle. For cosmic ray data, there are two indepen- 
dent trigger blocks: SciBar/EC and MRD. SciBar and 
the EC use a common cosmic ray trigger which is gen- 
erated using fast signals from the TA. The MRD has its 
own cosmic ray trigger which is also self-generated by 
discriminator outputs. Both SciBar/EC and the MRD 
collect 20 cosmic ray triggers in a cycle. 
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B. Detector simulation 



The GEANT4 framework is used for the detector simu- 
lation. The Bertini cascade model within GEANT4 
is used to simulate the interactions of hadronic parti- 
cles with detector materials. The detector simulation 
includes a detailed geometric model of the detector, in- 
cluding the detector frame and experimental hall and soil, 
which is based on survey measurements taken during de- 
tector construction. 

In the detector simulation of SciBar, low level data 
parameters are used as input to the simulation whenever 
possible. The energy loss of a charged particle in a single 
strip is simulated by GEANT, and this energy scale is 
tuned using cosmic ray data. Scintillator quenching is 
simulated using Birk's law [HOl with a value of Birk's 
constant, measured for K2K, of 0.0208 cm/MeV [H. 
The energy deposited by a charged particle is converted 
to photoelectrons using conversion factors measured for 
each channel with cosmic muons. The measured light 
attenuation length of each fiber (approximately 350 cm 
on average) is used in the simulation. Crosstalk between 
nearby MA-PMT channels is simulated using values mea- 
sured in a test stand prior to installation. The number 
of photoelectrons is smeared by Poisson statistics, and 
the single photoelectron resolution of the MA-PMT is 
simulated. To simulate the digitization of the PMT sig- 
nal, the number of photoelectrons is converted to ADC 
counts, and then electronics noise and threshold effects 
of the TA are simulated. 

TDC hit simulation includes light propagation delays 
in the WLS fibers. A logical OR of 32 MA-PMT channels 
is made for each TDC channel, and the time of each hit 
is converted to TDC counts. Multiple TDC hits in each 
channel are simulated. 

In the EC detector simulation, true energy deposition 
in scintillating fibers in the detector is converted to the 
number of photoelectrons using a conversion factor which 
is measured for each channel with cosmic-ray muons. The 
attenuation of light in the fiber is simulated using the 
measured attenuation length value. The number of pho- 
toelectrons is smeared by Poisson statistics and by the 
PMT resolution, and then converted to ADC counts. The 
time-dependent ADC gain due to the overshoot of the 
PMT signal is simulated based on a measurement with 
cosmic muons. Electronics noise is also simulated. 

For the detector simulation of the MRD, true energy 
deposition in each scintillator is converted to ADC counts 
using the conversion factor measured with cosmic muons. 
The attenuation of light in the scintillator as well as elec- 
tronics noise are simulated. Gaps between scintillator 
counters in each plane, which cause inefficiency, are in- 
cluded in the simulation. The time of energy deposition 
is digitized and converted into TDC counts. 



CM 
UJ 



0) 



o 
o 
a. 




Jun Jul Oct Nov Dec Jan Feb Mar Apr May Jun Jul Aug 

'"8 Date 




Oct Nov Dec Jan Feb Mar Apr May Jun Jul Aug 

Date 



FIG. 4: Experimental performance. In the top panel, the solid 
line shows the history of the accumulated number of protons 
on target and the dashed line shows the number of protons on 
target passing all data quality cuts. The bottom panel shows 
the number of charged current candidate events in SciBar 
normalized to the number of protons on target. The event 
rate difference between neutrino and antineutrino modes can 
be seen clearly. 



V. DATA SUMMARY 

The SciBooNE experiment took data from June 2007 
until August 2008. The data-taking is divided into three 
periods depending on the polarity of the horn, as sum- 
marized in Table HI Fig. IDJtop) shows a history of the 
accumulated number of protons on target; the two curves 
show the total protons on target for all events and the 
protons on target for events passing all data quality cuts, 
described below. 

In total, 2.64 x 10^° protons on target were delivered to 
the beryllium target during the SciBooNE data run. The 
beam datastream (measuring, for example, magnet cur- 
rent settings, measured beam intensity, measured peak 
horn current) is synchronized and merged with the corre- 
sponding SciBooNE detector datastream, provided that 
the spill time as measured by the beam instrumentation 
and by the detector match within 10 ms of each other. 
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TABLE I: Summary of SciBooNE data-taking. The table 
shows the number of protons on target (POT) collected after 
application of data quality cuts, as described in the text. 



Run Period 






POT 


Run 1 (Antineutrino) Jun. 2007 - 
Run 2 (Neutrino) Oct. 2007 - 
Run 3 (Antineutrino) Apr. 2008 - 


Aug. 
Apr. 
Aug. 


2007 
2008 
2008 


0.52 X 10^" 
0.99 X 10^° 
1.01 X 10^° 



Only spills that satisfy certain beam quality cuts are 
used for analysis. We require that the beam intensity 
is at least 0.1 x 10^^ protons per spill, that the agree- 
ment between the two toroid readouts along the beam- 
line is within 10%, that the absolute peak horn current 
is greater than 170 kA, and that the targeting efficiency 
is greater than 95%. Overall, beam quahty cuts reject 
less than 1% of the total number of protons on target 
accumulated during the run. A somewhat larger frac- 
tion of protons on target is rejected because of detector 
dead time, yielding about a 95% efficiency to satisfy all 
(beam plus detector) data quality cuts. After all beam 
and detector quality cuts, 2.52 x 10^" protons on target 
are usable for physics analyses. 

In this analysis, the full neutrino data sample is used, 
corresponding to 0.99 x 10^" protons on target satisfying 
all data quality cuts, collected between October 2007 and 
April 2008. During that time, all detector channels were 
operational on beam triggers. The experimental stability 
is demonstrated in Fig.lDJbottom), which shows the num- 
ber of charged current event candidates per protons on 
target. In this figure, the event reconstruction is a simple 
track finder which is used only for operations related 
studies, and not for the analysis described in this paper. 
The figure illustrates the event rate difference between 
neutrino mode and antineutrino mode running. 

The antineutrino data sample collected before and af- 
ter the neutrino data-taking period is not considered in 
this analysis. 

VI. EVENT RECONSTRUCTION AND 
ANALYSIS 

A. Track Reconstruction 

The first step of the event reconstruction is to search 
for two-dimensional tracks in each view of SciBar using 
a cellular automaton algorithm [s^ . For tracking, the 
hit threshold is set to two photoelectrons, correspond- 
ing to approximately 0.2 MeV. Three dimensional tracks 
are reconstructed by matching the timing and z-edges 
of the two dimensional tracks. The timing difference 
between two two dimensional tracks is required to be 
less than 50 nsec, and the z-edge difference must be less 
than 6.6 cm for upstream and downstream edges. Re- 
constructed tracks are required to have at least three- 



layer penetration, and therefore the minimum length of 
a reconstructed track is 8 cm in the beam direction. Ac- 
cording to the MC simulation, 96% of charged current 
interactions in SciBar are reconstructed to have at least 
one track. 

To identify charged current events, we look for events 
in which at least one reconstructed track in SciBar is 
matched with a track or hits in the MRD. Such a track 
is defined as a SciBar-MRD matched track. The most 
energetic SciBar-MRD matched track in any event is con- 
sidered as a muon candidate. For matching a MRD track 
to a SciBar track, the upstream edge of the MRD track 
is required to be on either one of the first two layers 
of the MRD. The transverse distance between the two 
tracks at the first layer of the MRD must be less than 
30 cm. The requirement on the difference between track 
angles with respect to the beam direction is given by 
|6'mrd - 6*861 < 6'max, whcrc 6lmax IS 'A functiou of 
the length of the MRD track, varying between 0.4 radian 
and 1.1 radians. For track reconstruction in the MRD, 
at least two hit layers in each view are needed, and thus 
this matching method is used for tracks which penetrate 
at least three steel plates. If no MRD track is found, we 
extrapolate the SciBar track to the MRD and search for 
nearby contiguous hits in the MRD identifying a short 
muon track. For matching MRD hits to a SciBar track, 
the MRD hit is required to be within a cone with an aper- 
ture of ±0.5 radian and a transverse offset within 10 cm 
of the extrapolated SciBar track at the upstream edge 
of the MRD. The timing difference between the SciBar 
track and the track or hits in the MRD is required to be 
within 100 nsec. The matching criteria impose a muon 
momentum threshold of 350 MeV/c. 

B. Particle Identification 

The SciBar detector has the capability to distinguish 
protons from muons and pions using dE/dx. The particle 
identification variable, Muon Confidence Level (MuCL) 
is calculated as follows. The confidence level at each 
plane is first defined as the fraction of events in the ex- 
pected dE/dx distribution of muons above the observed 
value, {dE / dx)ohs- The expected dE/dx distribution of 
muons is obtained by using cosmic-ray muons. Each 
plane's confidence level is combined to form a total con- 
fidence level, assuming the confidence level at each layer 
is independent. The MuCL is calculated as 

MuCL = Pxyi^l^ (1) 

i=0 

where n is the number of planes penetrated by the track, 
P = nr=i ^^ij CILi is the confidence level at the i-th 
plane. 

Fig. [5] shows the dE/dx distributions of muon and 
proton enriched samples. The predicted distributions 
of true muon and proton tracks are shown as hatched 
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FIG. 5: (Color online) dE/dx of muon enriched sample (left) 
and proton enriched sample (right). 




FIG. 6: (Color online) MuCL of muon enriched (left) and 
proton enriched (right) samples. 



histograms. To select muon candidates for this study, 
we first select SciBar-MRD matched tracks. According 
to the MC simulation, the sample is 94.7% pure muons 
with a small contamination of protons and charged pi- 
ous. For proton candidates, we select the second track in 
a charged current quasi-clastic (CC-QE) scattering en- 
riched sample made by cutting on a kinematic variable 
described later. The fraction of protons in the sample 
is 92.1%, estimated with the MC simulation. The con- 
tamination of charged pions and muons are estimated to 
be 5.5% and 1.6%, respectively. Proton candidates are 
clearly separated from muon candidates. 

The MuCL distributions for the muon enriched sam- 
ple and the proton enriched sample are shown in Fig. [B] 
Tracks with MuCL greater than 0.05 are considered 
muon-like (or pion-like) and the others are classified as 
proton-like. The probability of misidcntification is esti- 
mated to be 1.1% for muons and 12% for protons, aver- 
aged over track length in the muon and proton enriched 
samples. 



date on the most upstream layer of SciBar to eliminate 
incoming particles due to neutrino interactions in the up- 
stream wall or soil. The hit threshold for this veto cut is 
set to two photoelectrons. The neutrino interaction ver- 
tex is reconstructed as the upstream edge of the muon 
track. The vertex resolution is approximately 0.5 cm 
in each dimension, estimated with the MC simulation. 
We select events whose vertices are in the SciBar fidu- 
cial volume, defined to be ±130 cm in both the x and 
y dimensions, and 2.62 cm< z <157.2 cm, a total mass 
of 10.6 tons. The background contamination due to neu- 
trino events which occur in the EC and MRD is 2.0% and 
0.5%, respectively. Finally, the time of the muon candi- 
date is required to be within a 2 /isec window around 
the beam pulse. The cosmic-ray background contamina- 
tion in the beam timing window is only 0.5%, estimated 
using a bcam-ofi' timing window. According to the MC 
simulation, the selection efficiency and purity of true 
charged current events are 27.9% and 92.8%, respectively. 
Impurity comes from neutral current events (3.0%), 
77^ charged current events (1.6%), and neutrino events 
which occur in the EC/MRD (2.5%). The average neu- 
trino beam energy for true charged current events in the 
sample is 1.2 GeV. This SciBar-MRD matched sample 
is our standard charged current data set and defines the 
MC normalization, i.e. the MC distributions are normal- 
ized to the number of SciBar-MRD matched events in 
data. 

Two sub-samples of the SciBar-MRD matched sample 
are further defined; the MRD stopped sample and the 
MRD penetrated sample. Events with the muon stop- 
ping in the MRD are classified as MRD stopped events, in 
which we can measure the muon momentum. Events with 
the muon exiting from the downstream end of the MRD 
are defined as the MRD penetrated sample, in which we 
can measure only a part of the muon momentum. The 
average neutrino beam energy for true charged current 
events in the MRD stopped and MRD penetrated sam- 
ples are 1.0 GcV and 2.0 GeV, respectively, enabling a 
measurement of charged current coherent pion produc- 
tion at two different mean neutrino energies. 

The slopes of the muon angles with respect to the beam 
in the two SciBar views are used to calculate the three 
dimensional muon angle with respect to the beam {9^). 
The kinetic energy of the muon is calculated by the range 
and expected energy deposition per unit length (dE/dx) 
in SciBar, the EC and the MRD, 



kin 



^SB 

dE 
dx 



^EC 
^SB 



SB 



^MRD 

AEj^ 
cos 9,, 



E 



MRD 



MRD) 



(2) 



C. Charged Current Event Selection 

Events with at least one SciBar-MRD matched track 
are selected as charged current event candidates. We 
reject events with hits associated with the muon candi- 



where E^^, E^'~^, and E^^^ are the energy deposition 
in each detector. Lsb and Lmrd are the track length of 
the muon in SciBar and the range in the MRD, respec- 
tively. We set dE/dx\sB to 2.04 MeV/cm, and AE^'^, 
which is the energy deposited in the EC by a horizontally 
traversing minimum ionizing particle, is set to 91 MeV, 
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estimated with the GEANT4 simulation. £;MRD jg 
culated from a range to energy lookup table based on 
the MC simulation. For muons stopping in the MRD, 
the average muon momentum and muon angular reso- 
lutions are 50 MeV/c and 0.9 degree, respectively. For 
muons exiting the MRD, only a lower limit on muon mo- 
mentum is obtained, while the muon angle is determined 
with the same resolution as that of stopping muons. The 
systematic uncertainty in the muon momentum scale is 
estimated to be 2% which is dominated by the difference 
among various calculations of the range to energy lookup 
table. 



D. Event Classification 



• DATA 
^§ CC coherent it 
CC resonant n 
~ Other 
I I CCQE 




4 5 
Number of tracks 



The MRD stopped and MRD penetrated samples are 
further divided into sub-samples with the same selection 
criteria. Once the muon candidate and the neutrino in- 
teraction vertex are reconstructed, we search for other 
tracks originating from the vertex. For this purpose, 
the track edge distance is defined as the 3D distance 
between the vertex and the closer edge of another re- 
constructed track. Tracks whose edge distance is within 
10 cm are called vertex-matched tracks. Fig.[7]shows the 
distribution of the number of tracks at the vertex for the 
MRD stopped sample. For the MC simulation, the con- 
tributions from charged current coherent pion, charged 
current resonant pion, charged current quasi-elastic, and 
other interactions are shown separately. Most events are 
reconstructed as either one track or two track events. 
The two track sample is further divided based on parti- 
cle identification. We first require that the MuCL of the 
SciBar-MRD matched track is greater than 0.05 to reject 
events with a proton penetrating into the MRD. Then 
the second track in the event is classified as a pion-like 
or a proton-like track with the same MuCL threshold. 
Fig. [5] shows the contributions to the second track from 
true proton, pion, muon, and electron tracks as predicted 
by the MC simulation. 

In a charged current resonant pion event, up — > 
jj,~piT~^, the proton is often not reconstructed due to 
its low energy, and such an event is therefore identified 
as a two track ^ + tt event. To separate charged cur- 
rent coherent pion events from charged current resonant 
pion events, additional protons with momentum below 
the tracking threshold are instead detected by their large 
energy deposition around the vertex, so-called vertex ac- 
tivity. We search for the maximum deposited energy in 
a strip around the vertex, an area of 12.5 cm x 12.5 cm 
in both views. Fig. [9] shows the maximum energy for 
fi + TT events in the MRD stopped sample. A peak around 
6 MeV corresponds to the energy deposited in the strip 
containing the vertex by two minimum ionizing parti- 
cles, and a high energy tail is mainly due to the low 
energy proton. To simulate such protons, we consider 
re-interactions of nucleons in the nucleus as described 
in Sec. IIII El as well as ones outside the nucleus (de- 



FIG. 7: (Color online) Number of vertex-matched tracks for 
the MRD stopped sample. The MC distribution shown here 
is before tuning. 



Tc-like 




MuCL 



FIG. 8: (Color online) MuCL for the second track in the 
two-track sample. The MC distribution shown here is before 
tuning. 



scribed in Sec. lIVBp . De-excitation gamma-rays from 
the carbon nucleus do not affect the distribution since 
most of the gamma-rays first interact outside the ver- 
tex region. Events with energy deposition greater than 
10 MeV are considered to have activity at the vertex. 
Charged current coherent pion candidates are extracted 
from the sample oi fj, + tt events without vertex activ- 
ity. Four sub-samples, the one track events, fJ.+p events, 
fi+TT events with vertex activity, and fi+n events without 
vertex activity in the MRD stopped sample are used for 
constraining systematic uncertainties in the simulation, 
described next. 
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scale (i?pscaic), and a charged current quasi-elastic Pauli- 
suppression parameter k. All parameters are ratios to 
nominal values in the MC simulation, i.e. all parameters 
are set to 1 in the default MC simulation. 

The parameters i?2trfe/itrfci Rp/in and i?act represent 
possible migrations between subsamples due to system- 
atic uncertainties. The parameter i?pscaic changes the 
scale of the reconstructed muon momentum for the MC 
simulation. The parameter k, which was first intro- 
duced by MiniBooNE [s^l , controls the strength of Pauli- 
blocking and thus suppresses low charged current 
quasi-clastic events. We employ this parameter in the 
fitting because a deficit of data is found at low in the 
one track sample where the charged current quasi-clastic 
interaction is dominant. 

The function to be minimized is given by: 



FIG. 9: (Color online) Maximum deposited energy in a strip 
around the vertex for the ^ + vr events. The MC distribution 
shown here is before tuning. 



The term 
defined as 



iist 



X 



Xdist 



Asys 



(5) 



is calculated using a binned likelihood 



E. Tuning the Monte Carlo Simulation 

The MC simulation includes systematic uncertainties 
due to the detector response, nuclear effects, neutrino 
interaction models, and neutrino beam spectrum, and 
these uncertainties affect background estimation. The 
sources of systematic uncertainty are summarized in Sec- 
tion IVTlBj In order to constrain these uncertainties, the 
MC distributions of the square of the four-momentum 
transfer (Q^) arc fitted to the distributions of the four 
aforementioned data samples. The reconstructed is 
calculated as 



Ql, = 2El^^{E^ - cos e^)-ml 



(3) 



Xdist 



obs 



^obs 



(6) 



where P{n,v) = v'^e "/n! is the Poisson probability of 
finding n events with a expectation value v.^ and 



N^^^ are the observed and expected number of events 
in the i-th bin in subsamplc j (j =:onc track, ^ + p, 
jj, + TT with high and low vertex activity), respectively. 
The expected number of events for each sample is given 
by: 



where E'^°'^ is the reconstructed neutrino energy calcu- 
lated by assuming charged current quasi-clastic kinemat- 
ics, 



El, 



1 (rrip — m^) — (m 
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y)2 + 2S^(m„-l/) 



(m„ - V) - Ef,+Pf, cos 0f, 
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^cxp 



^norm ' J^2trk/ltrk 



'^i>7rH + R'rcs1ll°^^ii + -Rotlier«°*^^H 



^cxp 

I, fLTTL 



Rn 



-R2trk/ltrk ' 



where rup and m„ are the mass of proton and neutron, 
respectively, and V is the nuclear potential, which is set 
to 27 MeV. The one track events, n+p events, and ^ + tt 
events with and without vertex activity are fit simulta- 
neously. Each Q^.^c distribution is fit in bins of width 
0.05 (GcV/c)2 up to 1 (GeV/c)2. 

Wc introduce eight fitting parameters; the normaliza- 
tion factor of the MRD stopped sample (-Rnorm), the res- 
onant pion scale factor (i?res), the scale factor of other 
non-QE interactions (i?othcr), the ratio of the number 
of two track events to the number of one track events 
(i?2trk/itrk), the ratio of the number of ^ -I- p events to 

the number of (U-Ftt events (i?„/^), the ratio of the number -Rpscaie and k do not appear explicitly in these equa- 



QE 



i?resn'^°^ 



Mother ?^°^^^L 



(9) 



(10) 



are the number of charged cur- 



where n^^., nj^'^., 71°^^ 
rent quasi-clastic, charged current resonant pion, and 
other events in each bin in each subsamplc, respectively. 



of low vertex activity fi + ir events to the number of high 
vertex activity fi + n events (i?act), the muon momentum 



tions, but i?pscaic causes migration between Q bins and 
K changes n'^^j. 
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The term Xsvs' added to constrain systematic parame- 



ters, is calculated as: 

x'={Psys-Po)V-'{Psys 



(11) 



where Pgys represents the set of systematic parameters 
and Pq is the set of parameter values before fitting, ex- 
pressed as: 



sys 



R 



2trk/ltrk 



R, 



pscalc 



V is a covariancc matrix estimated by considering the 
possible variations due to systematic uncertainties in the 
detector responses, nuclear effects, neutrino interaction 
models, and neutrino beam spectrum. We prepare sev- 
eral MC event sets by changing each underlying physics 
parameter, i.e. the source of systematic uncertainty, by 
±1(7. The covariancc between two systematic parameters 
Pi and pj is calculated as: 



Vij = cov[p,,pj] = 



ApiApj\+ + ApiApjl^ 



(13) 



where ApiApj\_^(^_-^ is the product of variations of two 
parameters when the underlying physics parameter is in- 
creased (decreased) by the size of its uncertainty. The 
covariancc matrix is estimated to be: 



V = 



(0.20)2 
-(0.09)2 
f(0.10)2 




-(0.09)2 
(0.09)2 
-(0.07)2 




+ (0.10)2 
-(0.07)2 
(0.15)2 





(14) 



R 

norm ? Mother, -Ract, and K are unconstrained in the fit. 

Events with Q2^^ < g.lO (GeV/c)2 in the fi + n sam- 
ple with low activity are not included in the fit to avoid 
charged current coherent pion signal events. A data ex- 
cess is observed in the region with Q^^.^ < 0.15 (GeV/c)2 
in the fi+p sample. Further investigation reveals that the 
second track in the excess events is emitted at a relatively 
large angle with respect to the beam direction and has 
large dE/dx, thus the events have an additional large 
energy deposition at the vertex. Each of these events 
seems to have a muon and a proton with additional ac- 
tivity, and therefore the excess is not expected to affect 
the charged current coherent pion analysis. A possible 
candidate for the excess is charged current resonant pion 
production where the pion is absorbed in the nucleus. In 
such an event, two or more additional nuclcons should 
be emitted after the pion is absorbed, which is currently 
not simulated. The excess cannot be explained with the 
introduced fitting parameters, and therefore the region is 
not used in the fit. 

Fig. [To] shows reconstructed after the fitting for 
the one track, /i + p^ and /i + tt events with and with- 
out vertex activity. The best fit values and errors of the 
fit parameters are summarized in Table |TT1 These same 



fit parameters are also applied to the MRD penetrating 
sample. The xVd-O.f before the fit is 473/75 = 6.31. 
The xVd.o.f after the fit is 117/67 = 1.75. Even after 
fitting, the reduced is relatively large, which indicates 
that the introduced parameters are not sufficient in fully 
reproducing the data. To take into account the incom- 
pleteness of our simulation, we enlarge the errors on the 
fitting parameters by a factor of ^/)^Jd^ol. 



(12) TABLE II: Best fit values and errors of the fitting parameters 



Parameter 


Value 


Error 


-Rnorm 


1.103 


0.029 


-R2trk/ltrk 


0.865 


0.035 


Rp / n 


0.899 


0.038 


Ra.ct 


0.983 


0.055 


Rpscale 


1.033 


0.002 


Rtcs 


1.211 


0.133 


Rothei 


1.270 


0.148 


K 


1.019 


0.004 



F. Charged Current Coherent Pion Event Selection 

Charged current coherent pion candidates are ex- 
tracted from both the MRD stopped and MRD pene- 
trated samples with the same selection criteria. In this 
section, we describe the event selection for the MRD 
stopped sample. The event selection for the MRD pene- 
trated sample is summarized later. 

After selecting fi + tt events which do not have ver- 
tex activity, the sample still contains charged current 
quasi-elastic events in which a proton is misidentified 
as a minimum ionizing track. We reduce this charged 
current quasi-elastic background by making use of kine- 
matic information in the event. Since the charged cur- 
rent quasi-elastic interaction is a two-body interaction, 
one can predict the proton direction from the measured 
muon momentum p^ and muon angle 0^ ; 



Pp = i-Pf,x, -Pf,y, K""" - -Pm cos 0^ 



(15) 



where p^^x and are the projected muon momentum 
in the x and y dimension, respectively. E^J"^ is the recon- 
structed neutrino energy given by Equation |4l For each 
two-track event, we define an angle called A9p as the 
angle between the expected proton track direction given 
by Equation 1151 and the observed second track direction. 
Fig. [TT] shows the A6p distribution for /z + tt events in the 
MRD stopped sample. Events with A9p larger than 20 
degrees are selected. With this selection, 48% of charged 
current quasi-elastic events in the fi + tt sample are re- 
jected, while 91% of charged current coherent pion events 
pass the cut according to the MC simulation. 

Further selections are applied in order to separate 
charged current coherent pion events from charged cur- 



13 



-,5000 



> 

a 



a) 1 -track 



r4ooo 



U3000 



2000 



1000 




DATA 

CC coherent 
CC resonant tt 
Other 
CC QE 



, 500 



> 

a 



(b) i^+p 



400 



.5! 300 



200 



100 




• DATA 

g CC coherent tt 
^ CC resonant tt 

Other 
■ CCOE 



■a^ 500 




c) p,+7c with activity 

• DATA 

% CC coherent tt 
CC resonant tt 
Other 
CC QE 



DATA 

CC coherent it 
CC resonant n 
Other 




20 40 60 80 100 120 140 160 180 

Ae (degrees) 



FIG. 11: (Color online) AOp for the /i + 7r events in the MRD 
stopped sample after fitting. 
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FIG. 12: (Color online) Track angle of the pion candidate 
with respect to the beam direction for the + tt events after 
the charged current quasi-elastic rejection after fitting. 




(GeV/cf 

FIG. 10: (Color online) Reconstructed after fitting for (a) 
the one track, (b) IJ.+P, (c) jj. + -K with activity, and (d) /i + vr 
without activity samples. 



rent resonant pion events which are the dominant back- 
grounds for this analysis. Fig. [T2l shows the angular dis- 
tribution of pion candidates with respect to the beam 
direction. In the case of charged current coherent pion 
events, both the muon and pion tracks are directed for- 
ward. Events in which the track angle of the pion can- 
didate with respect to the beam direction is less than 
90 degrees arc selected. 

Fig. [T2] shows the reconstructed distribution for 
the fj, + n events after the pion track direction cut. Al- 
though a charged current quasi-elastic interaction is as- 
sumed, the of charged current coherent pion events 
is reconstructed with a resolution of 0.016 (GeV/c)^ and 
a shift of -0.024 (GeV/c)^ according to the MC simu- 
lation. Finally, events with reconstructed less than 
0.1 (GeV/c)^ are selected. The charged current coherent 
pion event selection is summarized in Table lllll In the 



14 




Q' (GeV/c)' Q' (GeV/c)' 



FIG. 13: (Color online) Reconstructed for the /^ + 7r events 
in the MRD stopped sample after the pion track direction cut 
and after fitting. 



FIG. 14: (Color online) Reconstructed for the /^ + 7r events 
in the MRD penetrated sample after the pion track direction 
cut after fitting. 



TABLE III: Event selection summary for the MRD stopped 
charged current coherent pion sample. 



Event selection 


DATA 


MC 


Coherent tt 






Signal 


B.G. 


Efficiency 


Generated in SciBar fid.vol. 




1,939 


156,766 


100% 


SciBar-MRD matched 


30,337 


978 


29,359 


50.4% 


MRD stopped 


21,762 


715 


20,437 


36.9% 


2 track 


5,939 


358 


6,073 


18.5% 


Particle ID + tt) 


2,255 


292 


2,336 


15.1% 


Vertex activity cut 


887 


264 


961 


13.6% 


CC-QE rejection 


682 


241 


709 


12.4% 


Pion track direction cut 


425 


233 


451 


12.0% 


Reconstructed cut 


247 


201 


228 


10.4% 



signal region, 247 charged current coherent pion candi- 
dates are observed, while the expected number of back- 
ground events is 228±12. The error comes from the errors 
on the fitting parameters summarized in Tabic [TTl The 
background in the final sample is dominated by charged 
current resonant pion production. The "other" back- 
ground is comprised of 50% charged current DIS, 32% 
neutral current, and 18% 77^ events. The selection effi- 
ciency for the signal is estimated to be 10.4%. 



G. MRD penetrated Charged Current Coherent 
Pion Events 



The same selection is applied to the MRD penetrated 
sample to extract charged current coherent pion candi- 
dates at higher energy. Fig. [14] shows the reconstructed 
distribution of the MRD penetrated charged cur- 
rent coherent pion sample. The reconstructed and 
for the MRD penetrated sample are calculated from 



TABLE IV: Event selection summary of MRD penetrated 
charged current coherent pion sample. 



Event selection 



DATA MC 

Signal B.G. 



Coherent tt 
Efficiency 



Generated in SciBar fid.vol. 
SciBar-MRD matched 



MRD penetrated 
2 track 

Particle ID (/x -I- n) 
Vertex activity cut 
CC-QE rejection 
Pion track direction cut 
Reconstructed cut 





1,939 156,766 


100% 


30,337 


978 


29,359 


50.4% 


3,712 


177 


4,375 


9.1% 


1,029 


92 


1,304 


4.7% 


418 


78 


474 


4.0% 


167 


71 


186 


3.6% 


134 


67 


135 


3.5% 


107 


66 


109 


3.4% 


57 


60 


40 


3.1% 



muon angle and partially-reconstructed muon energy, us- 
ing Equation [3] and Equation [H respectively. Although 
only a part of the muon energy is observed, the re- 
construction performance is essentially same because of 
the small muon angle. The event selection is summarized 
in Table ITVl In the signal region, 57 charged current co- 
herent pion candidates arc observed, while the expected 
number of background events is 40±2.2. The background 
in the final sample is dominated by charged current res- 
onant pion production. The "other" background is com- 
prised of 75% charged current DIS, and 25% events. 
The selection efficiency for the signal is estimated to be 
3.1%. 
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VII. RESULTS 



Cross Section Ratio 



1. MRD stopped charged current coherent pion sample 

After subtracting background and correcting for the 
selection efficiency, tlie number of cliarged current co- 
herent pion candidates in the MRD-stopped sample is 
measured to be 179 ± 190(stat); this error includes the 
uncertainty in the background estimation. No evidence 
of charged current coherent pion production is found in 
the sample. The neutrino energy dependence of the se- 
lection efficiency for charged current coherent pion events 
is shown in Fig. 1151 The mean neutrino beam energy for 
true charged current coherent pion events in the sample 
is estimated to be 1.1 GeV after accounting for the ef- 
fects of the selection efficiency. The RMS of the neutrino 
beam energy is 0.27 GeV. 

The total number of charged current interactions is 
estimated by using the SciBar-MRD matched sample. 
We observe 30,337 SciBar-MRD matched events. As de- 
scribed in section IVIi the selection efficiency and purity 
of charged current events are estimated to be 27.9% and 
92.8%, respectively. The neutrino energy dependence 
of the selection efficiency for charged current events is 
shown in Fig. [161 After correcting for the efficiency and 
purity, the number of charged current events is measured 
to be (1.091 ± 0.006(stat)) x 10^. 

Using this information, the ratio of the charged current 
coherent pion to total charged current production cross 
sections is measured to be (0.16 ± 0.17(stat)to'27(sys)) x 
10~^ at 1.1 GeV, where the systematic error is described 
later. The result is consistent with the non-existence 
of charged current coherent pion production, and hence 
we set an upper limit on the cross section ratio by us- 
ing the likelihood distribution (C) which is convolved 
with the systematic error. We calculate the 90% con- 
fidence level (C.L.) upper hmit (UL) using the relation 

J^^ Cdx/ Cdx = 0.9 to be: 

a(CC coherent vr) ^ ^ ^ 
(t(CC) 

at a mean neutrino energy of 1.1 GeV. 



(16) 



fects of the selection efficiency. The RMS of the neutrino 
beam energy is 0.80 GeV. 

Due to the higher neutrino energy in the charged cur- 
rent coherent pion sample, the MRD penetrated charged 
current sample is chosen to estimate the number of to- 
tal charged current interactions at a similar neutrino en- 
ergy. We observe 3,712 MRD penetrated events, and 
the efficiency and purity of true charged current 
events are estimated to be 4.5% and 97.5%, respec- 
tively. The impurity largely comes from charged cur- 
rent events. After correcting the efficiency and purity, 
the number of charged current events is measured to be 
(0.804 ± 0.013(stat)) x 10^. A 26% difference between 
the MRD matched and penetrated samples is found while 
the estimated uncertainty due to the neutrino flux is 14%. 
However, this is expected to be a small effect on the cross 
section ratio measurement. 

The ratio of the charged current coherent pion to total 
charged current production cross sections is measured to 
be (0.68±0.32(stat)ti^:^^(sys)) x 10"^ at 2.2 GeV. The 
systematic error is described later. No significant ev- 
idence for charged current coherent pion production is 
observed, and hence we set an upper limit on the cross 
section ratio at 90% C.L.: 

^(CC coherent tt) , 9 

a(CC) < ^ (1^) 

at a mean neutrino energy of 2.2 GeV. 



B. Systematic Uncertainties 

The sources of systematic error are divided into five 
categories, (i) detector response and track reconstruction, 
(ii) nuclear effects, (iii) neutrino interaction models, (iv) 
neutrino beam, and (v) event selection. We vary these 
sources within their uncertainties and take the resulting 
change in the cross section ratio as the systematic un- 
certainty of the measurement. Table IVl summarizes the 
uncertainties in the charged current coherent pion cross 
section ratio for the MRD stopped and MRD penetrated 
samples. The total systematic error is 2? ^ 10 for 
the MRD stopped sample, and +0 25 x 10"^ foj. ^he MRD 
penetrated sample. 



Detector Response and Track Reconstruction 



2. MRD penetrated charged current coherent pion sample 

After subtracting background and correcting for the se- 
lection efficiency, the number of charged current coherent 
pion candidates in the MRD penetrating sample is mea- 
sured to be 548 ± 254(stat). As in the MRD stopping 
sample, this includes the uncertainty due to the back- 
ground estimation. The mean neutrino beam energy for 
true charged current coherent pion events in the sample 
is estimated to be 2.2 GeV after accounting for the ef- 



The crosstalk of the MA-PMT was measured to be 
3.15% for adjacent channels, with an absolute error of 
0.4%. The single photoelectron resolution of the MA- 
PMT is set to 50% in the simulation, to reproduce the 
observed dE/dx distribution of cosmic muons. The ab- 
solute error is estimated to be ±20%. Birk's constant 
of the SciBar scintillator was measured to be 0.0208 ± 
0.0023 cm/MeV [5l| and is varied within the measure- 
ment error to evaluate the systematic. The hit threshold 
for track reconstruction is varied by ±20%. 
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CC coherent % 




TABLE V: Summary of the systematic errors in the charged 
current coherent pion cross section ratio. 



Source MRD stopped 

error (xlO~^) 


MRD penetrated 
error (xl0~^) 


Detector response +0.10 —0.18 
Nuclear effect +0.20 -0.07 
Neutrino interaction model +0.17 —0.04 
Neutrino beam +0.07 -0.11 
Event selection +0.07 -0.14 


+0.18 -0.18 
+0.19 -0.09 
+0.08 -0.04 
+0.27 -0.13 
+0.06 -0.05 


Total +0.30 -0.27 


+0.39 -0.25 



MRD stopped 

+ 



MRD penetrated 



Ev (GeV) 

FIG. 15: (Color online) Neutrino energy spectra and selection 
efhciencies as a function of neutrino energy for charged current 
coherent pion events. 




Ev (GeV) 
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FIG. 16: (Color online) Neutrino energy spectra and selection 
efficiencies as a function of neutrino energy for all charged 



current events. 



2. Nuclear Effects 

We consider uncertainties in final state interactions in- 
side the nucleus. This includes rescattering of nucleons 
and pions in the initial target nucleus. For pions pro- 
duced by neutrino interactions, uncertainties on the cross 
sections for pion absorption and pion inelastic scattering 
in the nucleus are considered. The cross section of pion 
charge exchange is negligible compared with the other ef- 
fects and is hence neglected. In the momentum range of 
pions from A decays, the cross section measurement un- 
certainty for both absorption and inelastic scattering is 
approximately 30% [i^l ■ Nucleon re-scattering in the nu- 
cleus affects vertex activity. The uncertainty in the cross 
section is estimated to be 10%. In the NEUT simulation, 
the Fermi momentum of nucleons is set to 217 MeV/c 
for carbon. According to electron quasi-elastic scatter- 



ing data [5JI, the value is approximately 221 ± 5 MeV/c. 
Therefore, an uncertainty of ±5 MeV/c is assigned. 



3. Neutrino Interaction Models 

In the NEUT simulation, we set the axial vector mass 
Ma to 1.21 GeV/c^ for both QE and resonant pion pro- 
duction. The uncertainty in this value is estimated to be 
approxim ately ±0.1 GeV/c^ based on recent measure- 
ments [3lll32|: results from past experiments are system- 
atically lower than the recent measurements [55|, and 
thus we only vary Ma to 1.11 GeV/c^, and take that 
change as the systematic error. We consider the uncer- 
tainty in the charged current resonant /x~n7r+//i~p7r+ 
cross section ratio. The uncertainty in this ratio is esti- 
mated to be 7% using four SciBooNE data samples de- 
scribed in Section IVIDI In addition, a disagreement of 
the shape is observed in the /i + tt events with vertex 
activity where charged current resonant pion production 
is dominant, as shown in Fig. [101 We estimate the sys- 
tematic uncertainty, due to evident low suppression 
of charged current resonant pion production even after 
the MC tuning, by re-weighting the true distribution 
of charged current resonant pion events according to the 
observed low deficit. 
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4- Neutrino Beam 

The uncertainties in secondary particle production 
cross sections in proton-beryllium (p-Be) interactions, 
hadronic interactions in the target or horn, and the horn 
magnetic field model are varied within their externally 
estimated error bands. Detailed descriptions of each un- 
certainty are found elsewhere [23|. Uncertainties asso- 
ciated with the delivery of the primary proton beam to 
the beryllium target and the primary beam optics, which 
result in the overall normalization uncertainty, are not 
considered in this analysis since it cancels in the cross 
section ratio. 



larger in antineutrino data because of the reduced to- 
tal 77 charged current event rate. Because of this, the 
antineutrino search has the potential to be even more 
sensitive. 

Furthermore, theoretical models also make concrete 
connections between the charged and neutral current co- 
herent pion production processes. As mentioned in Sec- 
tion m the MiniBooNE Collaboration has already pub- 
lished an observation of neutral current coherent pion 
production in the same neutrino beam as SciBooNE. The 
SciBooNE neutral current coherent pion search is, there- 
fore, also interesting and may shed considerable light on 
the behavior of this interaction process. 



5. Event Selection 

For the event selection variable A^p, we evaluate the 
systematic uncertainty in the cross section ratio by vary- 
ing the cut placement. The uncertainty in the A9p cut 
for charged current quasi-elastic rejection is estimated 
to be ±5 degrees. For the other variables, we already 
tuned the MC simulation using the migration parameters 
or considered possible systematic sources. Therefore, we 
do not include additional systematic uncertainties due to 
these selections. 



C. Discussion 

Having not observed evidence for charged current co- 
herent pion production, we set 90% C.L. upper limits 
at two different mean neutrino energies. According to 
the Rein-Sehgal model 0, Q implemented in our simu- 
lation, the cross section ratio of charged current coher- 
ent pion production to total charged current interactions 
is expected to be 2.04 xlO~^. Our limits correspond 
to 33% and 67% of the prediction at 1.1 GeV and 2.2 
GeV, respectively. For reference, the total charged cur- 
rent cross section averaged over the MRD stopped and 
MRD penetrated samples are 1.05 x 10~^® cm^/nucleon 
and 1.76 x 10"'^* cm^/nucleon, respectively, estimated 
with the MC simulation. Our results are consistent with 
the K2K result; cr(CC coherent 7r)/cr(CC) < 0.60 x 10"^ 
at 90% C.L. measured in a 1.3 GeV wide-band neutrino 
beam. As shown in Fig. [2l several recent models pre- 
dict a considerably smaller coherent cross section, which 
appears consistent with our results. 

Because of the connection of the neutrino and antineu- 
trino coherent pion production processes in the theoret- 
ical models, it will be interesting to repeat this analysis 
on SciBooNE 's already collected antineutrino data. Most 
models predict a similar absolute cross section for neu- 
trino and antineutrino coherent pion production, which 
means the ratio of charged current coherent pion events 
to charged current inclusive events is expected to be 



VIII. CONCLUSIONS 

In conclusion, we have searched for muon neutrino 
charged current coherent pion production on carbon in 
the few GeV region using the full SciBooNE neutrino 
data set of 0.99 x 10^" protons on target. No evidence 
of charged current coherent pion production is found, 
and hence we set 90% C.L. upper limits on the cross 
section ratio of charged current coherent pion to total 
charged current production cross sections at 0.67 x 10"^ 
and 1.36 x 10~^, at mean neutrino energies of 1.1 GeV 
and 2.2 GeV, respectively. 
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